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The  operation  of  a  dual-free-layer  (DFL)  spin-torque  nano-oscillator  (STNO)  is  studied  is  demon¬ 
strated  that  in  a  practically  interesting  regime  when  the  magnetizations  of  the  two  free  layers  (FLs) 
precess  in  opposite  directions  along  large-angle  out-of-plane  trajectories,  thus  doubling  the  gener¬ 
ation  frequency,  the  operation  of  the  DFL  STNO  is  strongly  hysteretic  as  a  function  of  a  bias  dc 
current.  The  stable  magnetization  dynamics  starts  at  a  rather  large  magnitude  of  the  bias  dc  current 
density  Jdc  >  j(][sh  when  the  bias  current  is  increased,  but  the  regime  of  stable  counter-precession 
of  the  FLs  persists  till  rather  low  magnitudes  of  the  bias  dc  current  density  Jt’j(w  <  Jdc  <  j(][sh  when 
the  bias  current  is  decreased.  This  hysteresis  is  caused  by  the  dipolar  coupling  between  the  FLs, 
and  is  especially  pronounced  for  small  distances  between  the  FLs  and  small  magnetic  damping  in 
them.  The  discovered  hysteretic  behavior  of  the  DFL  STNO  implies  the  possibility  of  application 
of  a  strong  initial  pulse  of  the  bias  current  (greater  than  the  upper  threshold  j(][gh  of  the  stable 
dynamics)  and  subsequent  reduction  of  the  bias  current  to  a  working  point  (J*j(w  <  Jdc  <  J^gh) 
corresponding  to  the  required  output  frequency  /(Jdc).  The  obtained  results  might  be  important 
for  the  practical  development  of  DFL  STNOs  with  optimized  operation  characteristics. 

PACS  numbers:  85.75.-d,  07.57.Hm,  84.40.-x 


I.  INTRODUCTION 

The  spin-transfer  torque  (STT)  effect  [1-4]  in  mag¬ 
netic  multilayers  theoretically  predicted  by  L.  Berger  and 
J.  C.  Slonczewski  provides  a  new  method  of  manipulation 
of  the  magnetization  dynamics  in  nano-magnetic  systems 
[5,  6].  The  spin-torque  nano-oscillators  (STNOs),  based 
on  the  STT  effect,  are  now  considered  as  promising  base 
elements  for  various  nano-scale  microwave  devices,  par¬ 
ticularly,  microwave  signal  sources  [7-17]  and  microwave 
detectors  [18-25].  However,  for  these  devices  to  become 
practical  in  modern  microwave  signal  processing  systems 
it  is  necessary  to  enlarge  the  range  of  their  operational 
frequencies,  to  improve  their  phase  noise  characteristics, 
and  to  increase  their  integrated  output  power  [26,  27].  It 
is,  also,  preferable  that  these  devices  should  be  able  to 
operate  in  a  zero  (or  very  low)  bias  dc  magnetic  field. 

A  step  towards  the  development  of  such  a  practical 
STNO  has  been  made  in  [28],  where  an  oscillator  with  a 
perpendicular  polarizer  has  been  suggested  and  demon¬ 
strated  experimentally.  In  the  STNO  proposed  in  [28] 
the  free  magnetic  layer  (FL)  was  situated  between  the 
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perpendicular  polarizer  and  an  in-plane  magnetized  ref¬ 
erence  layer.  The  magnetization  of  the  free  layer  was 
pushed  out  of  plane  by  STT  from  the  polarizer  and  at 
a  sufficiently  large  bias  current  it  was  precessing  along 
a  large  amplitude  trajectory  around  the  resulting  easy- 
plane  shape  anisotropy  field  [29,  30].  The  stationary  pre¬ 
cession  of  the  FL  magnetization  with  respect  to  the  static 
perpendicular  magnetization  of  the  reference  layer  results 
in  the  microwave-frequency  oscillations  of  the  device  re¬ 
sistance  transformed  in  a  microwave  voltage  signal  by 
means  of  the  giant  magneto-resistance  (GMR)  effect  be¬ 
tween  the  layers  [28].  Because  of  the  relatively  large  os¬ 
cillation  amplitude  this  type  of  STNO  offers  a  significant 
improvement  in  microwave  power  emission  over  conven¬ 
tional  STNOs  consisting  of  an  in-plane  FL  and  an  in¬ 
plane  pinned  magnetic  layer  (PL).  The  main  limitation 
of  this  structure,  however,  is  a  relatively  low  maximum 
frequency  caused  by  the  formation  of  a  static  vortex  in 
the  FL  at  high  densities  of  the  dc  current  [29] . 

In  this  paper,  we  consider  a  symmetric  dual-free-layer 
(DFL)  STNO,  containing  two  perpendicularly  polarized 
PLs  and  two  adjacent  FLs  with  easy-plane  anisotropy 
(Fig.  1(a)).  It  has  been  suggested  recently  that  such 
a  DFL  STNO  could  be  used  as  a  high-frequency  mi¬ 
crowave  source  in  zero  bias  magnetic  field,  because  in 
such  a  geometry  it  is  possible  to  achieve  the  opposite  di- 
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FIG.  1:  (Color  Online)  (a)  Layout  of  the  symmetric  dual-free- 
layer  (DFL)  STNO  consisting  of  four  magnetic  layers:  two 
inner  free  magnetic  layers  (FLs)  of  the  thickness  L  and  two 
outer  pinned  layers  with  perpendicular  polarization  separated 
by  the  GMR  spacer  of  the  thickness  d.  (b)  The  torques  acting 
on  the  magnetization  vectors  of  the  FLs  in  the  DFL  STNO 
driven  by  the  transverse  bias  dc  current  of  the  density  Jdc- 
The  conservative  torque  (Tp)  caused  by  the  effective  mag¬ 
netic  field  Beff  induces  magnetization  precession  in  the  FLs 
around  the  normal  to  the  STNO  layers.  The  non-conservative 
spin-torques  (Ts)  caused  by  the  perpendicular  polarizers  ad¬ 
jacent  to  each  of  the  FLs  tend  to  increase  the  precession  angle, 
while  the  non-conservative  torques  (Tpj)  caused  to  the  Gilbert 
damping  in  the  FLs  tend  to  reduce  it. 


rections  of  the  current-driven  magnetization  precession 
in  the  FLs,  and,  therefore,  to  double  the  generated  mi¬ 
crowave  frequency  (the  angle  between  the  magnetizations 
of  the  FLs  in  a  DFL  DTNO  changes  twice  as  fast  as  the 
angle  between  the  magnetizations  of  the  free  and  pinned 
layers)  [31].  Below  we  consider  the  regime  of  a  DFL 
STNO  operation,  when  the  FL  magnetizations  precess 
in  the  opposite  directions  along  the  large-angle  out-of- 
plane  (OOP)  trajectories  and  generate  a  relatively  large 
microwave  power  at  a  frequency  that  is  a  sum  of  the 
magnetization  precession  frequencies  in  individual  FLs. 

Although  the  physical  mechanism  of  the  frequency 
doubling  in  the  proposed  regime  of  the  DFL  STNO  op¬ 
eration  seems  to  be  simple  and  obvious,  it  is  not  unique. 
It  has  been  shown  in  [32]  that  the  frequency  doubling  in 
a  DFL  STNO  structure  can  be  also  achieved  when  the 
magnetizations  of  the  FLs  precess  in  the  same  direction, 
but  are  phase-shifted.  In  our  theoretical  considerations 
below  we  will  not  explicitly  consider  this  “phase-shifted” 
case.  Recent  micromagnetic  simulations  performed  in 
[33]  have  shown  that  the  regime  of  the  DFL  STNO  oper¬ 
ation  with  counter-precessing  FLs  is  limited  by  a  rather 
high  dc  current  density  threshold  Jt h  needed  for  the  exci¬ 
tation  of  a  stable  OOP  precession  in  the  FLs  ( Jdc  >  Jt h) 
and  by  the  formation  of  a  vortex  in  the  FL  at  higher 
dc  current  densities  (Jdc  >  Jvortex)-  It  is  important  to 
note,  that  the  dc  current  density  threshold  Jth  found 
from  the  micromagnetic  simulations  is  of  the  order  of 
Jth  ~  0.35  •  10s  A/cm2  for  a  DFL  STNO  with  typical 
parameters  [33]  and  substationally  exceeds  the  threshold 
for  a  conventional  single- free-layer  STNO  with  similar  pa¬ 


rameters.  Consequently,  the  range  of  the  DFL  STNO  op¬ 
eration  frequencies  in  the  regime  with  counter-precessing 
FLs  (~  Jvortex/ Jth)  is  smaller  than  the  similar  range  for 
a  conventional  single- free-layer  STNO.  This  limited  fre¬ 
quency  range  complicates  the  practical  applications  of 
the  DFL  STNO  for  microwave  signal  generation,  while 
a  rather  high  dc  current  density  threshold  Jth  limits  the 
reliability  of  such  devices. 

In  this  work,  using  a  simple  numerical  macrospin 
model  described  in  Section  II,  we  study  the  operation 
of  a  DFL  STNO  in  the  regime  with  counter-precessing 
FLs  in  a  wide  range  of  bias  dc  current  densities  Jdc  and 
reveal  the  existence  of  a  strong  hysteresis  in  the  depen¬ 
dence  of  the  operational  parameters  of  a  DFL  STNO  on 
the  bias  current  density.  In  the  Section  III  we  show  that 
a  stable  counter-precession  of  magnetizations  of  the  FLs, 
indeed,  starts  at  a  rather  large  magnitude  of  the  bias 
current  density  Jdc  =  Jt*[[sh  (approximately  equal  to  the 
bias  dc  current  density  threshold  Jt h  found  from  the  mi¬ 
cromagnetic  simulations  [33] )  when  the  bias  current  den¬ 
sity  is  increased.  However,  the  regime  of  stable  counter¬ 
precession  persists  till  rather  low  magnitudes  of  the  bias 
current  Jdc  =  Jt1£w  when  the  bias  current  density  is  de¬ 
creased.  We  demonstrate  that  the  lower  threshold  Jt*]]w 
is  mainly  determined  by  the  damping  in  the  system  (like 
the  threshold  in  a  conventional  single-free-layer  STNO), 
while  the  upper  threshold  Jt*]]sh  of  the  auto-oscillation 
regime  is  mainly  determined  by  the  existence  of  a  dipo¬ 
lar  coupling  between  the  FLs  and/or  by  the  in-plane 
anisotropy  in  the  FLs. 

Therefore,  a  relatively  large  external  bias  current  is 
needed  only  to  start  the  auto- oscillations  in  the  device, 
and  after  the  auto-oscillations  have  started  it  is  possible 
to  significantly  reduce  the  bias  current  density  to  Jt']]w 
while  preserving  the  regime  of  stable  counter-precessing 
magnetizations  of  the  FLs  of  the  device  at  a  working 
point  (./]/"  <  Jdc  <  Jm®11)  corresponding  to  the  required 
output  frequency  /(Jdc)-  We  believe  that  the  obtained 
results  might  be  important  for  the  practical  development 
of  DFL  STNOs  with  optimized  working  characteristics. 


II.  MODEL 

A.  Schematic  and  fundamentals  of  the  DFL  STNO 

The  layout  of  a  symmetric  DFL  STNO  is  shown  in 
Fig.  1(a).  The  STNO  is  a  circular  nano-pillar  consisting 
of  four  metallic  ferromagnetic  layers  separated  by  three 
non-magnetic  spacers.  The  outermost  magnetic  layers 
serve  as  current  polarizers,  and  are  composed  of  the  ma¬ 
terials  with  perpendicular  magnetic  anisotropy  that  is 
sufficiently  large  to  keep  the  magnetizations  of  these  lay¬ 
ers  out  of  the  sample  plane  along  the  z-axis  of  the  co¬ 
ordinate  system  shown  in  Fig.  1(b)  (p*  =  z  are  the  unit 
vectors  in  the  direction  of  the  magnetizations  of  the  po¬ 
larizers,  i  =  t  and  i  =  b  for  the  top  and  bottom  layers, 
respectively,  z  is  the  unit  vector  of  the  z-axis). 
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The  magnetizations  Mj  of  the  two  inner  ferromag¬ 
netic  layers  lie  in  the  plane  of  the  sample  at  equilib¬ 
rium,  and  these  layers  serve  as  the  FLs  of  the  STNO 
(Mj  =  Mfti  m i ,  Mfti  is  the  saturation  magnetization  of 
the  zth  FL,  m;  is  the  unit  vector  along  the  magnetization 
vector  in  the  zth  FL,  i  =  {t,  b}). 

In  the  presence  of  a  direct  bias  current  applied  across 
the  entire  nano-pillar,  as  it  is  shown  in  Fig.  1(b),  the  per¬ 
pendicular  polarizers  transfer  the  spin  torque  to  the  adja¬ 
cent  FLs,  Ts,i  ~  m,  x  (mj  x  p.j),  and  this  torque  pushes 
the  magnetizations  of  the  FLs  out  of  the  (x,y)  plane 
towards  the  opposite  ends  of  the  structure  (m2jb  <  0, 
mZjt  >  0;  ?n2j  j  is  ^-component  of  mj).  The  torque 
Td,i  ~  (m,  x  (dnii/dt))  created  by  Gilbert  damping 
opposes  this  action  of  the  spin  torque  and  pushes  the  FL 
magnetizations  back  towards  the  plane  of  the  device.  The 
easy-plane  demagnetizing  fields  arise  in  both  FLs  in  pro¬ 
portion  to  their  out  of  plane  magnetization  components 
(Bd)i  ~  — m2|»z),  and,  as  shown  in  Fig.  1(b),  induce  the 
FLs’  magnetizations  to  undergo  precession  in  opposite 
directions  around  the  normal  to  the  sample  plane  [34]  in 
a  zero  bias  dc  magnetic  held.  This  precession  due  to  the 
GMR  effect  produces  oscillations  of  the  resistance  of  the 
whole  nano-pillar,  and  the  frequency  of  these  oscillations 
is  equal  to  the  sum  of  the  individual  precession  frequen¬ 
cies  of  the  FLs.  Thus,  using  the  above  described  design 
of  a  DFL  STNO  it  is  possible  to  approximately  double 
the  generation  frequency  attainable  with  a  single  FL  and 
fixed  reference  layer.  The  amplitude  of  the  resistance  os¬ 
cillations  in  a  DFL  STNO  is  similar  to  the  full  amplitude 
of  the  GMR  between  the  FLs,  resulting  in  a  relatively 
large  microwave  power  emitted  from  the  device  [33]. 


B.  A  simplified  model  of  a  DFL  STNO 

Although  the  layout  of  a  DFL  STNO  presented  in  Sec¬ 
tion  II A  contains  four  magnetic  layers  (two  inner  FLs 
and  two  outer  PLs  (see  Fig.  1(a))  in  the  following  we  as¬ 
sume  that  the  magnetizations  of  the  PLs  of  the  STNO 
are  completely  fixed  and  lie  strictly  perpendicular  to  the 
layers’  plane,  and,  therefore,  we  consider  only  the  STT- 
induced  dynamics  in  the  two  FLs.  We  use  the  macrospin 
approximation  to  describe  the  magnetization  dynamics 
in  the  FLs,  i.e.  we  assume  that  each  of  the  dimensionless 
vectors  irq  =  M,  /M ^ ,  directed  along  the  magnetiza¬ 
tion  of  the  zth  FL  Mj  depends  only  on  time  t  for  both 
FLs  (z  =  {t,  b}).  We  also  assume  that  the  considered 
system  operates  in  a  zero  bias  dc  magnetic  field  and  is 
completely  symmetric  (the  thicknesses  and  magnetic  pa¬ 
rameters  of  both  FLs  are  the  same.  We  consider  a  case 
when  the  FLs  of  the  STNO  are  coupled  with  one  another 
by  a  dipole-dipole  interaction  and  are  characterized  by 
an  easy-axis  in-plane  anisotropy,  the  magnitude  of  which 
will  vary  in  our  simulations. 

Despite  the  obvious  limitations  of  the  chosen  model  it 
is  sufficiently  accurate  in  the  case  when  the  density  Jdc  of 
the  bias  dc  current  is  lower  ( J<jc  <  Jvortex)  than  the  den¬ 


sity  Jvortex  corresponding  to  the  vortex  formation  thresh¬ 
old  calculated  in  [33].  We  show  below  in  Section  III  that 
our  simplified  numerical  macrospin  model  of  the  DFL 
STNO  allows  one  to  correctly  predict  the  value  of  the 
threshold  Jt h  corresponding  to  the  appearance  of  stable 
dynamics  determined  from  the  much  more  complicated 
and  time-consuming  micromagnetic  simulations  [33].  It 
also  satisfactory  describes  the  dependence  of  the  gener¬ 
ated  microwave  frequency  on  the  density  of  the  bias  dc 
current  Jdc  for  current  densities  below  the  threshold  of 
vortex  formation  [33]. 


C.  The  equations  for  the  magnetization  dynamics 
in  a  DFL  STNO 


The  dynamics  of  the  normalized  magnetization  vectors 
rrij  (t)  =  Mj  ( t)/Mf  in  the  zth  FL  (z  =  {t,  b})  under  the 
action  of  a  bias  dc  current  with  density  Jdc  is  governed  by 
the  Landau-Lifshitz-Gilbert  equation  with  an  additional 
Slonczewski-Berger  STT  term  describing  the  influence  of 
a  spin  transfer  torque  [1,  2]: 

drrij 

—JT  =  +  TD,i  +  T s,i  •  (1) 

at 

The  first  term,  T  pj ,  in  the  right-hand  side  part  of 
this  equation  describes  the  torque  created  by  the  effec¬ 
tive  magnetic  field  B0ff;j  that  causes  the  conservative 
precession  of  the  unit  vector  mj  about  the  direction  of 
this  effective  magnetic  field  [35]: 

TP,j  =7[Beff)i  x  mj]  .  (2) 

Here  7  =  gy, b/H  ss  2tt  ■  28  GHz/T  is  the  modulus  of  the 
gyromagnetic  ratio,  g  is  the  spectroscopic  Lande  factor, 
/.zb  is  the  Bohr  magneton,  and  h  is  the  reduced  Planck 
constant. 

The  second  term,  is  the  dissipative  Gilbert 

torque  describing  the  energy  dissipation,  which  can  be 
written  as  [35]: 


T 


D,i 


=  a 


mj  x 


drrij 

dt 


(3) 


where  a  is  the  dimensionless  Gilbert  damping  parameter 
of  the  FL  (the  same  for  both  FLs). 

The  last  term  in  the  right-hand  side  part  of  Eq.  (1), 
Tg>  j ,  is  the  Slonczewski-Berger  STT  term  describing  the 
interaction  of  the  magnetization  with  the  spin-polarized 
current  traversing  zth  FL.  The  expression  for  this  torque 
obtained  in  [1,  2]  can  be  written  as: 


Ts,j  =  —  <Tj  Jdc,  i  [mj  x  (mj  x  p> )]  ,  (4) 

where  Jdc,  j  is  the  bias  dc  current  density  traversing  the 
STNO  (connected  to  the  bias  dc  current  /cic,  i  by  the  ex¬ 
pression  /dc,j  =  7rf?2Jdc,ij  where  R  is  the  radius  of  the 
nano-pillar),  pj  is  the  unit  vector  along  the  magnetiza¬ 
tion  of  the  PL  (pj  =  z  in  our  simulations),  and  <jj  is  the 
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“current  density-spin  torque”  proportionality  coefficient, 
which  can  be  written  in  a  general  form  as: 

cq  =  c>o?(P,  A) .  (5) 


Here  op  =  (•yK/2e)/(MfL)  is  the  angular-independent 
factor  in  the  “current  density-spin  torque”  coefficient  cq  , 
e  is  the  modulus  of  the  electron  charge,  Mf  is  the  satu¬ 
ration  magnetization  of  the  FL,  A  is  the  thickness  of  the 
FL,  and  ?(P,  A)  is  the  factor  describing  the  angular  de¬ 
pendence  of  STT  depending  on  the  dimensionless  degree 
P  (0  <  P  <  1)  of  the  current  spin-polarization  and  the 
angle  A  between  the  magnetizations  of  ith  FL  and  its  cor¬ 
responding  polarizer.  The  angular  dependence  c(P,  A)  is 
different  for  the  different  spacer  materials  between  the 
magnetic  layers.  In  the  case  of  a  metallic  spacer  (con¬ 
sidered  in  this  work)  this  factor  can  be  approximated  by 
[1]: 


<r(AA) 


'(1  +  P)3 

4p3/2 


(3  +  cos  A) 


(6) 


where  cos  A  =  (m*  •  Pi)-  We  specify  the  current  density 
JdCj  j  in  Eq.  (4)  as: 

t  _  /  +  Jdc  ,  i  =  t  (Top  free  layer)  ,  , 

dc’*  (  — Jdc  ,  i  =  b  (Bottom  free  layer)  ’  '  ' 


In  Eqs.  (10)  we  took  into  account  the  fact  that  the  FLs 
have  identical  parameters,  and,  therefore,  the  dimension¬ 
less  tensor  elements  k  and  p  depending  on  the  geometry 
of  the  structure  are  the  same  for  both  the  top  and  bottom 
FLs. 

In  the  practically  interesting  case  of  thin  magnetic  lay¬ 
ers,  when  the  layers’  thickness  L  and  the  thickness  of  the 
GMR  spacer  between  the  FLs  d  are  substantially  smaller 
than  the  radius  of  the  structure  H(L,d<i?),  the  tensor 
elements  n  and  p  are  given  by  the  following  expressions 
[36]: 


K 

P 


1  A 
2nR 


C  —  In 


C_L 
2 nR 

-f( 


1 

~~  4tt  (L/R) 


—  F 


d  +  L 
R 


(11a) 


(lib) 


where  C  =  ln(8)  —  1/2  and  F{ £)  =  (2  In  |£|.  Typically,  n 
and  p  are  small  quantities  ( k ,  p  <C  1).  For  instance,  in  the 
case  of  the  DFL  STNO  with  chosen  typical  parameters 
(see  Section  II D)  k  =  0.0707,  p  =  0.01134. 

An  easy-axis  anisotropy  field  of  the  ith  FL  was  chosen 
to  be  directed  along  the  x-axis,  and  is  defined  as  [35] 


with  opposite  signs  in  the  two  FLs  due  to  the  different 
order  in  which  the  spin-polarized  electrons  traverse  the 
two  polarizer/free-layer  pairs. 

In  our  model  the  effective  magnetic  field  Beff. 
in  Eq.  (2)  has  the  contributions  from  the  self¬ 
demagnetization  held  B *  j,  cross-demagnetization  held 
B,  j  (the  sum  of  these  helds  is  the  total  demagnetiza¬ 
tion  held  acting  on  ith  FL;  i,  j={t,  b}),  and  an  easy-axis 
anisotropy  held  Ba,i  : 

Beff,  i  =  Bj,  i  +  B ij  +  Ba,  i  ■  (8) 

It  is  important  to  note,  that  the  effective  magnetic  held 
Bcff,  i  does  not  include  an  external  bias  dc  magnetic  held, 
which  stresses  one  more  time  the  possibility  of  normal 
operation  of  a  DFL  STNO  in  a  zero  applied  dc  held. 

The  self-demagnetization  held  B,;,j  and  the  cross¬ 
demagnetization  held  B  ij  are  dehned  as  follows  [36]: 

B,  i  =  — /x0M/Nj,  jiiij ,  (9a) 

B  ,  =  ~Po M/Njjnij,  (9b) 

where  /i0  is  the  vacuum  magnetic  permeability,  N,  is 
the  self-demagnetization  tensor 

(  k  0  0  \ 

Nj,j  =  0  k  0  ,  (10a) 

\  0  0  1-2 k  J 

and  Njj  is  the  cross-demagnetization  tensor: 

(  P  0  0  \ 

N,,j  =  I  0  p  0  .  (10b) 

\  0  0  -2p  J 


Ba,  i  —  Ba  (mj  •  x)  x ,  (12) 

where  Ha  >  0  is  the  magnitude  of  the  anisotropy  held 
along  the  easy  x-axis  in  both  FLs. 

The  details  about  the  numerical  procedure  of  our 
macrospin  simulations  are  presented  in  the  Appendix  (see 
Section  V  A  for  details). 

D.  Parameters  of  the  DFL  STNO 

For  a  quantitative  description  of  the  DFL  STNO’s  op¬ 
eration  we  used  the  following  typical  parameters  of  the 
nano-pillar  structure  [33].  We  considered  a  symmetric 
circular  nano-pillar  with  the  radius  R  =  25  nm,  3  nrn 
thick  FLs  (L  =  At  =  Ab  =  3  nm)  and  12  nm  thick  GMR 
spacers  ( d  =  12  nm)  (see  Fig.  1(a)). 

We  assumed  that  both  FLs  have  the  saturation  mag¬ 
netization  Mf  =  Mf  t  =  Mf  b  =  560  emu/ern3  (thus 
po Mf  «  704  mT),  which  is  typical  for  thin  Permalloy 
(Nii-^FeA  hlms  sandwiched  between  the  Cu  layers  [37]. 

We  assumed  a  value  of  the  current  polarization  to  be 
P  =  0.3,  which  is  typical  for  GMR  spin  valves.  It  has 
been  shown  recently  in  [33]  that  the  value  of  the  Gilbert 
damping  parameter,  a,  has  a  strong  impact  on  the  dy¬ 
namics  of  a  DFL  STNO,  and,  therefore,  we  simulated  the 
dynamic  behavior  of  our  model  system  in  a  wide  range 
of  a  values  0.002-0.09  (the  same  for  both  FLs).  The  sub¬ 
stantial  increase  of  the  FL  damping  can  be  realized  in 
practice  by  doping  the  FL  with  Tb  [38].  We  also  found 
that  the  magnetization  dynamics  in  a  DFL  STNO  signif¬ 
icantly  depends  on  the  magnitude  of  the  magnetic  field 
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FIG.  2:  (Color  Online)  (a)  Numerically  calculated  depen¬ 
dences  of  the  precession  frequencies  in  the  top  mx,t  (ft)  and 
bottom  mx,b  (/b)  FLs  in  a  DFL  STNO  on  the  density  Jdc  of 
the  bias  dc  current,  demonstrating  strong  hysteresis.  Dynam¬ 
ics  starts  at  Jdc  =  J^gh  when  the  bias  current  increases  (red 
squares  and  cyan  circles)  and  stops  at  Jdc  =  J*hW  when  the 
current  decreases  (orange  hollow  squares  and  blue  crosses). 
The  system  can  be  in  the  in-plane  static  (IPS)  or  in  the  out- 
of-plane  precessional  (OPP)  state,  (b)  Numerically  calculated 
dependence  of  the  frequency  of  an  output  signal  /  in  a  DFL 
STNO  on  the  bias  direct  current  density  Jac.  The  gener¬ 
ated  frequency  is  equal  to  the  sum  of  the  precession  frequen¬ 
cies  in  the  FLs,  f  =  ft  +  fb  ■  Yellow  star  marks  a  possible 
working  point,  which  could  be  reached  by  using  the  proce¬ 
dure  described  in  Section  IIIC.  The  curves  are  calculated  for 
a  =  0.08,  Ba  =  5  mT.  All  the  other  parameters  of  the  DFL 
STNO  are  presented  in  Section  II D. 


Ba  of  the  easy-axis  in-plane  anisotropy  (see  Eq.  (12)), 
and,  therefore,  performed  simulations  in  which  the  field 
Ba  (again,  the  same  for  both  FLs)  was  varied  in  the 
range  0-20  mT. 


III.  RESULTS  AND  DISCUSSION 

A.  Hysteresis  in  the  dependence  of  the  generated 
frequency  on  the  bias  current 


Using  the  macrospin  simulation  procedure,  described 
in  details  in  Section  V  A,  we  numerically  solved 
the  coupled  Landau-Lifshitz-Gilbert  equations  with 
Slonczewski-Berger  STT  terms  (1)  for  the  magnetizations 
of  the  two  FLs  and  calculated  the  frequencies  of  magneti¬ 
zation  precession  in  the  top  and  bottom  FLs.  The  typical 
dependence  of  the  frequency  generated  in  a  DFL  STNO 


on  the  density  Jdc  of  the  applied  bias  dc  current  is  shown 
in  Fig.  2.  As  one  can  see,  the  frequency  dependence  on 
the  bias  current  is  different  depending  on  whether  the 
current  density  increases  or  decreases,  i.e.  this  depen¬ 
dence  demonstrates  a  strong  hysteresis. 

When  the  density  of  the  dc  bias  current  is  increased 
(red  squares  and  cyan  circles  in  Fig.  2(a),  red  squares 
in  Fig.  2(b)),  the  magnetizations  in  both  FLs  stay  in 
plane  (IPS)  without  any  precession  till  the  current  den¬ 
sity  reaches  a  rather  large  magnitude  J^sh  ( Jt^gh  ss 
0. 33-10®  A/cm2  in  Fig.  2).  Above  this  critical  (threshold) 
value  the  magnetizations  in  both  FLs  start  to  precess  in 
opposite  directions,  and  the  frequency  of  precession  in¬ 
creases  with  the  increase  of  the  current  density.  Note, 
that  the  generated  frequency  just  above  the  threshold  has 
a  rather  large  value  (about  8  GHz,  see  Fig.  2(b)).  How¬ 
ever,  when  the  density  of  the  bias  dc  current  is  decreased 
(orange  hollow  squares  and  blue  crosses  in  Fig.  2(a),  blue 
hollow  points  in  Fig.  2(b)),  the  regime  of  stable  counter¬ 
precession  of  the  FLs  persists  till  rather  low  magnitudes 
of  the  bias  current  density  Jt'£w  (Jjgw  «  0.12- 10®  A/cm2, 
(see  Fig.  2),  and  the  generated  frequency  which  decreases 
with  the  decrease  of  current  can  be  much  lower  (down  to 
approximately  2  GHz). 

Note,  that  the  threshold  value  J^8*1  (at  which  the  sta¬ 
ble  counter-precession  in  two  FLs  starts)  obtained  in  our 
macrospin  simulations  is  very  close  to  the  threshold  value 
Jth  ~  0.35  •  10®  A/cm2  determined  from  the  micromag- 
netic  simulations  [33].  This  is  a  clear  confirmation  that 
our  simplified  “macrospin”  model  is  sufficiently  accurate 
to  get  even  the  quantitative  information  about  the  char¬ 
acteristics  of  DFL  STNO.  One  can  also  note,  that  the 
difference  between  the  frequencies  of  magnetization  pre¬ 
cession  in  the  top  and  bottom  layer  of  the  DFL  STNO 
(ft  and  /b,  respectively)  is  caused  by  the  angular  depen¬ 
dence  of  the  STT  given  by  Eq.  (6).  If  in  our  simulations 
we  assume  that  cos/?*  =  const  in  Eq.  (6),  the  hysteresis 
behavior  of  the  DFL  STNO  remains  the  same  as  before, 
but  the  precession  frequencies  of  the  top  and  bottom  lay¬ 
ers  ft  and  fb  become  equal. 

The  origin  of  the  hysteresis  in  DFL  STNO  may  be  qual¬ 
itatively  explained  as  follows.  In  the  absence  of  a  bias 
dc  current  the  system  exists  in  the  so-called  “in-plane 
state”  (IPS).  In  this  state  the  magnetizations  of  the  FLs 
are  dipolarly  coupled  to  each  other  and  the  energy  of  the 
system  is  minimal  when  the  magnetizations  have  oppo¬ 
site  directions.  In  other  words,  the  dipolar  coupling  cre¬ 
ates  a  potential  well  for  the  magnetizations  (the  in-plane 
anisotropy  has  a  similar  effect),  and  a  finite  torque  (and, 
therefore,  a  finite  bias  current  density  J^1®11)  is  required 
to  move  the  magnetizations  out  of  this  potential  well. 
Thus,  the  higher  threshold  current  is  determined 

by  the  depth  of  this  potential  well  and  depends,  mainly, 
on  the  strengths  of  the  dipolar  coupling  between  the  FLs 
and/or  on  the  strength  of  the  in-plane  anisotropy  field. 
However,  once  the  dipolar  coupling  between  the  two  FLs 
is  destroyed  by  a  sufficiently  strong  spin  torque  created 
by  the  bias  current,  a  much  smaller  spin  torque  is  suffi- 
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FIG.  3:  (Color  Online)  Numerically  calculated  dependences  of  the  current  density  thresholds  (blue  circles  and  line)  and 
^thSh  (red  squares  and  line)  on  the  different  parameters  of  the  DFL  STNO.  (a)  The  dependence  of  j’™  and  Jt^sh  on  the 
Gilbert  damping  parameter  a.  (b)  The  dependence  of  J*£w  and  Jt^gh  on  the  magnitude  of  the  in-plane  anisotropy  field  Ba- 
(c)  The  dependence  of  j'£w  and  J^sh  on  the  distance  between  the  FLs  d.  Green  stars  and  dashed  line  show  the  dependence  of 
the  cross-demagnetization  coefficient  p  on  d,  p(d).  All  the  curves  are  calculated  for  a  =  0.08,  Ba  =  5  mT  (if  other  values  of 
parameters  are  not  mentioned  specifically).  All  the  other  parameters  of  the  DFL  STNO  are  presented  in  Section  II D. 


cient  to  sustain  the  persistent  magnetization  precession: 
as  in  the  other  types  of  STNO,  the  spin  torque  only  needs 
to  compensate  the  natural  Gilbert  damping  in  the  FLs. 
Thus,  the  lower  critical  current  Jt^w  is  determined  by  the 
energy  losses  and  should  be  approximately  proportional 
to  the  Gilbert  damping  constant  a. 


B.  Threshold  bias  current  densities  of  the  DFL 
STNO 

To  verify  the  above  described  qualitative  picture,  we 
performed  detailed  investigations  of  the  dependence  of 
the  critical  current  densities  Jt^sh  and  j'£w  on  various 
system  parameters.  In  accordance  with  the  above  pre¬ 
sented  explanations,  the  critical  currents  in  DFL  STNO 
should  depend  on  three  major  factors:  (a)  Gilbert  damp¬ 
ing  parameter  a  of  the  system,  (b)  easy-axis  anisotropy 
field  with  magnitude  Ba,  and  (c)  dipolar  coupling  be¬ 
tween  the  FLs,  dependent  on  the  system  geometric  pa¬ 
rameters  (e.g  the  distance  d  between  the  FLs)  and  char¬ 
acterized  by  the  cross-demagnetization  coefficient  p(d). 
The  numerically  calculated  dependences  of  both  critical 
currents  on  these  parameters  are  presented  in  Fig.  3. 

One  can  see  from  Fig.  3(a),  that  the  upper  thresh¬ 
old  (red  squares  and  line)  does  not  depend  on  the 

damping  parameter  a,  while  the  lower  threshold  ,/t’£w 
(blue  circles  with  line)  increases  linearly  with  a.  The 
dependence  Jt^w(a)  is  similar  to  a  similar  relation  for 
a  conventional  single-free-layer  STNO,  where  the  mag¬ 
netization  of  a  FL  can  precess  relative  “freely” ,  when  an 
applied  bias  current  compensates  the  natural  damping  in 
the  system.  In  that  case  in  a  conventional  STNO  there 
is  a  linear  relation  between  the  threshold  current  and  the 
Gilbert  damping  parameter  a  and  a  similar  dependence 
JyJw(a)  is  seen  in  our  DFL  STNO. 


At  the  same  time,  the  fact  that  the  higher  threshold 
j‘;;sh  is  independent  of  damping  shows  that  the  origin 
of  this  upper  threshold  current  is  non-dissipative.  In 
contrast  the  curve  of  a  linear  dependence  of  on 

the  magnitude  of  the  in-plane  anisotropy  field  Ba  (red 
squares  and  line),  shown  in  Fig.  3(b),  indicates  that  a 
strong  in-plane  anisotropy  field  acts  in  the  same  man¬ 
ner  as  a  dipolar  coupling  between  the  FLs:  it  prevents 
the  appearance  of  a  stable  magnetization  precession  and 
creates  a  finite  dc  current  density  threshold  for  this  pre¬ 
cession  which  increases  linearly  with  the  increase  of  the 
anisotropy  field.  In  contrast  to  the  clear  linear  depen¬ 
dence  Jt1^gh(I?A),  the  dependence  of  the  low  current  den¬ 
sity  threshold  Jjjjw  on  Ba  (blue  circles  and  line,  see 
Fig.  3(b))  is  relatively  weak. 

The  other  interesting  result  is  that  the  above  de¬ 
scribed  hysteresis  exists  even  in  the  absence  of  an  in¬ 
plane  anisotropy  field.  In  that  case  its  existence  is  deter¬ 
mined  by  the  dipolar  coupling  between  the  FLs  only.  The 
curves  Jt1(1lgh(d)  (red  squares  and  line),  Jt1£w(d)  (blue  cir¬ 
cles  and  line),  and  p{d)  (green  stars  and  line)  calculated 
for  different  values  of  the  distance  between  the  FLs  d  are 
presented  in  Fig.  3(c)  These  curves  convincingly  show  a 
substantial  effect  of  the  dipolar  coupling  between  the  FLs 
on  the  magnetization  dynamics  in  the  DFL  STNO. 

Comparing  the  curves  Jt^sh(d)  and  p(d)  it  becomes 
obvious,  that  the  upper  current  density  threshold  is  pro¬ 
portional  to  the  cross-demagnetization  coefficient  of  the 
system  p,  Jt^gh  ~  p.  This  means  that  the  upper  current 
density  threshold  J^lgh  can  be  selected  by  choosing  an  ap¬ 
propriate  geometry  of  the  system  giving  some  particular 
value  of  the  parameter  p.  Although  the  low  current  den¬ 
sity  threshold  j\™  also  depends  on  p,  this  dependence  is 
substantially  weaker  in  comparison  to  the  dependence  of 
the  upper  current  density  threshold.  Consequently,  one 
can  easily  change  the  width  of  the  hysteresis  loop  defined 
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as  (,/tl||gh  —  Jt'hW)  by  choosing  a  particular  value  of  the  ge¬ 
ometric  factor  p  dependent  on  the  thickness  L  and  radius 
R  of  the  FLs  and  on  the  distance  d  between  the  FLs  (see 
Eq.  (lib)  for  details). 


C.  Stable  magnetization  dynamics  in  a  DFL  STNO 
for  bias  current  densities  below  the  upper  threshold 

The  discovered  hysteretic  behavior  of  the  DFL  STNO 
implies  the  possibility  of  application  of  a  strong  initial 
pulse  of  the  bias  current  (greater  than  the  upper  thresh¬ 
old  of  the  stable  dynamics  J^]8*1)  and  subsequent  re¬ 
duction  of  the  bias  current  to  a  working  point  corre¬ 
sponding  to  the  required  output  frequency.  In  order  to 
prove  this  statement  we  performed  a  numerical  experi¬ 
ment  when  the  STNO  is  biased  by  a  dc  current  density 

JdcW  =  jrk + 

The  supplied  bias  current  signal  consists  of  two  inde¬ 
pendent  component:  a  constant  “working”  dc  current 
with  density  J™°rk  =  0.15  •  10s  A/cm2,  that  is  substan¬ 
tially  smaller  than  the  upper  threshold  of  the  stable  dy¬ 
namics  J^j]8'1  (the  working  point  corresponding  to  the 
J™rk  is  shown  by  the  yellow  star  in  Fig.  2),  and  an  ini¬ 
tial  pulse  of  the  bias  current  characterized  by  a  current 
density  J^“lse(t)  and  duration  Tpuise. 

In  the  numerical  experiment  we  studied  how  the  mag¬ 
netizations  of  the  FLs,  characterized  by  the  unit  vec¬ 
tors  mt(f),  mb(t),  evolve  from  the  initial  equilibrium 
anti-parallel  configuration  [mt(0)  =  (1,0,0),  nib(0)  = 
(—1,0,0)]  over  time  t,  when  the  bias  current  density 
Jdc{t)  passes  through  the  DFL  STNO.  We  considered 
three  different  cases:  when  the  bias  current  density  is 
constant  (J^lse  =  0,  Fig.  4(a)),  when  the  rectangu¬ 
lar  pulse  of  the  bias  current  has  magnitude  of  J^lse  = 
0.1  •  108  A/cm2  not  exceeding  the  higher  threshold  of  the 
magnetization  dynamics  (Fig.  4(b)),  and  the  third  case 
when  jPculse  =  0.2  •  10s  A/cm2  is  over  the  higher  thresh¬ 
old  of  the  magnetization  dynamics  (Fig.  4(c)).  The  pulse 
duration  in  the  last  two  cases  was  T’puise=  0.5  ns. 

It  can  be  seen  from  Fig.  4  that  the  dynamics  of  the  sys¬ 
tem,  described  by  the  time  profiles  of  the  in-plane  mXt t  (t) 
(blue  curve)  and  out-of-plane  mz> t  (t)  (red  curve)  compo¬ 
nents  of  the  unit  vector  mt  (t)  in  the  top  FL  [39]  of  the 
DFL  STNO,  critically  depends  on  the  amplitude  J^lso 
of  the  pulsed  component  of  the  bias  current.  A  weak 
constant  dc  current  density  J™rk  induces  only  small  per¬ 
turbations  of  the  IPS  stationary  state  of  the  DFL  STNO 
(see  Fig.  4(a)).  When  the  initial  pulse  is  present  in  the 
bias  current,  but  the  total  initial  amplitude  of  the  cur¬ 
rent  density  ,/™rk  +  JP“lse  does  not  exceed  the  higher 
threshold  of  the  magnetization  dynamics  in  the  DFL 
STNO,  such  a  pulse  creates  a  transient  dynamical  pro¬ 
cess,  but  the  IPS  state  of  the  system  remains  stable  after 
the  end  of  the  current  pulse  (see  Fig.  4(b)).  The  situa¬ 
tion  changes  dramatically  when  a  strong  initial  pulse  of 
the  bias  {J^lse  =  0.2  •  108  A/cm2)  is  applied.  If  the  to¬ 


tal  initial  amplitude  of  the  current  density  J™rk  +  J^lse 
exceeds  the  higher  threshold  J^8’1  of  the  stable  magneti¬ 
zation  dynamics  in  the  DFL  STNO  the  stable  magnetiza¬ 
tion  precession  is  excited  in  both  FLs  of  the  DFL  STNO, 
and  this  dynamics  persists  after  the  end  of  the  initial  cur¬ 
rent  pulse  (see  Fig.  4(c)  for  the  magnetization  dynamics 
in  the  top  FL). 

The  obtained  results  clearly  show  that  it  is  possible 
to  achieve  a  regime  of  stable  counter-precession  of  FLs  in 
the  DFL  STNO  for  relatively  low  bias  dc  current  densities 
(below  the  the  higher  threshold  of  the  DFL  STNO  stable 
dynamic)  if  a  short  pulse  of  bias  current  exceeding  this 
threshold  is  applied  initially. 


IV.  SUMMARY 

In  conclusion,  we  report  that  in  a  DFL  STNO  oper¬ 
ating  in  the  absence  of  a  bias  dc  magnetic  field  there  is 
a  strong  hysteresis  in  the  dependence  of  the  frequency 
generated  by  the  STNO  on  the  bias  dc  current.  When 
the  bias  current  is  increased  the  microwave  dynamics  in 
the  DFL  STNO  starts  at  a  relatively  high  value  of  a  bias 
current  density,  but  when  the  bias  current  is  decreased 
this  dynamics  persists  till  a  much  lower  value  of  the  bias 
current  density.  The  lower  threshold  current  is  deter¬ 
mined  by  the  Gilbert  damping  parameter  (precessional 
state  vanishes  due  to  the  damping),  while  the  higher 
threshold  in  the  bias  current  is  mainly  caused  by  the 
existence  of  the  in-plane  anisotropy  and  dipolar  coupling 
between  the  FLs.  The  discovered  hysteretic  behavior  of 
a  DFL  STNO  implies  the  possibility  of  application  of  a 
sufficiently  strong  initial  pulse  of  the  bias  current  (larger 
than  the  higher  threshold  of  the  stable  magnetization  dy¬ 
namics)  and  subsequent  reduction  of  the  bias  current  to 
a  working  point  corresponding  to  the  required  generated 
frequency.  We  believe  the  obtained  results  will  be  impor¬ 
tant  for  the  practical  development  of  DFL  STNOs  with 
optimized  working  characteristics. 
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FIG.  4:  (Color  Online)  Numerically  calculated  time  profiles  of  the  in-plane  ( mX}t  ( t ),  blue  curve)  and  out-of-plane  (mZ} t  ( t ),  red 
curve)  magnetization  components  in  the  top  FL  of  the  DFL  STNO  for  different  regimes  of  time-dependent  current  biasing  (b,  d, 
f).  The  DFL  STNO  is  biased  by  a  dc  current  density  having  a  constant  and  a  pulsed  components  Jdc(t)  =  Jd°lk  +  and 

its  time  profile  is  shown  by  the  green  curve  with  green  shading  (a,  c,  e).  Dashed  horizontal  lines  show  the  levels  corresponding 
to  the  higher  (orange  lines)  and  lower  (violet  lines)  current  thresholds  J^w).  (a)  The  constant  current  density  below 

the  higher  threshold  JJ(nlk  =  0.15  •  108  A/cm2  induces  only  small  perturbations  of  the  IPS  magnetization  state  (b).  (c)  The 
current  density  containing  a  constant  part  and  a  pulse  of  the  duration  0.5  ns  with  the  amplitude  Jd/lse  =  0.1  •  108  A/cm2  (still 
below  the  higher  threshold)  induces  some  transient  dynamics,  but  in  the  end  the  IPS  magnetization  state  remains  stable  (d). 
(e)  The  current  density  containing  a  constant  part  (above  the  lower  threshold)  and  a  strong  pulse  of  the  duration  T,pUiae=0.5  ns 
with  the  amplitude  J^“lse  =  0.2  •  10s  A/cm2  (exceeding  the  higher  threshold)  induces  a  stable  magnetization  precession  (f). 
The  curves  are  calculated  for  a  =  0.08,  Ba  =  5  mT.  All  other  parameters  of  the  DFL  STNO  are  presented  in  Section  II D. 


V.  APPENDIX 
A.  Simulation  procedure 

Our  macrospin  simulations  were  based  on  the  numeri¬ 
cal  solution  of  a  pair  of  coupled  Landau-Lifshits-Gilbert- 
Slonczewski  (LLGS)  equations  (1)  for  both  top  and  bot¬ 
tom  FLs  biased  by  a  dc  current  density  (7),  where  the 
effective  magnetic  field  is  given  by  Eq.  (8)  and  the  spin 
torque  term  is  given  by  Eq.  (5). 

In  order  to  explore  the  spectrum  of  a  microwave  signal 
generated  by  the  DFL  STNO  under  the  action  of  a  bias 
dc  current,  we  assumed  the  current  density  c  to  be  a 
slow  variable  Jdc(t),  that  changes  with  time  linearly  at  a 


sufficiently  slow  rate: 

nt,  o  <t<^ 

,  ^<t<Tsim  ’ 
(13) 

where  77  =  \dJdc{t)/dt\  is  the  time  rate  of  the  current 
density  variation,  J™cax  =  0.57yTsim  is  the  maximum  den¬ 
sity  of  the  bias  dc  current  traversing  the  DFL  STNO 
and  Tsim  is  the  total  simulation  time.  During  the  first 
half  of  that  time,  0  <  t  <  0.5Tsjm ,  the  current  den¬ 
sity  increases,  while  during  the  second  half  of  that  time, 
0.5Tsim  <  t  <  Tsim,  it  decreases.  We  performed  a  series 
of  simulations  with  different  values  of  77  and  found  that 
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the  chosen  value  of  77,  rj  =  5- 107  A/ (cm2 /is),  does  not  no¬ 
ticeably  perturb  the  magnetization  dynamics,  so  that  the 
variations  of  the  current  density  (13)  can  be  considered 
as  quasi-static. 

While  solving  the  coupled  LLGS  equations  (1)  we  used 
the  following  initial  conditions  for  the  unit  vectors  at 
t  =  0:  mt(0)  =  (1,0,0),  mb(0)  =  (—1,0,0).  As  a  result 
of  a  numerical  solution  we  obtained  the  time  profiles  of 
the  unit  vectors  mt  =  mt(i)  and  mb  =  mb(t)  in  the 
range  of  time  0  <  t  <  Ts;m  in  top  and  bottom  FLs, 
respectively.  Using  the  discrete  Fourier  transform  pro¬ 
cedure  for  the  ^-components  of  the  calculated  profiles, 
mXjt  (t)  and  mXib(t),  we  calculated  the  frequencies  of  the 
magnetization  oscillations  in  the  top  and  bottom  FLs, 
ft  and  /b ,  respectively.  The  characteristic  accuracy  of 
the  calculated  frequencies  in  our  simulations  was  about 
50  MHz  ,  which  was  sufficient  for  the  quantitative  com¬ 


parison  of  our  results  with  the  results  of  micronragnetic 
simulations  and  possible  laboratory  experiments. 

When  analyzing  the  possibility  to  excite  a  stable  mag¬ 
netization  dynamics  in  the  FLs  for  bias  dc  current  densi¬ 
ties  Jdc  lower  than  the  upper  threshold  Jt^|s  ,  we  assumed 
that  the  DFL  STNO  is  biased  by  the  dc  current  density 
consisting  of  two  components.  The  first  one  was  an  ini¬ 
tial  current  pulse  of  the  density  magnitude  and 

duration  Tpuise,  while  the  second  one  was  a  constant  cur¬ 
rent  density  J™rk  which  determined  the  required  output 
frequency  of  a  microwave  signal  generated  by  the  DFL 
STNO: 


Jdc(t) 


Twork  Tpulse 
Jdc  f  'Me  J 


Twork 
Jdc  ? 


0  <  t  <  TpUise 

Tpulse  f  '  1  A  Tsim 


(14) 
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